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Abstract 

In spite of good performance of Tomlinson-Harashima 
Precoding (THP) in MIMO systems, it suffers from some 
drawbacks, including modulo, shaping and power loss. 
Moreover, in time-varying MIMO channels, due to channel 
variations, there is some more loss in performance of THP. 
To overcome this problem, the Channel State Information 
(CSI) should be updated each (or some) intervals. Sending 
channel information for each time slot needs a high volume 
of feedback data from receiver to transmitter. In temporal 
correlated MIMO channels, it is possible to use the 
correlation property of the channel and hence compensating 
channel uncertainty by using lower feedback volume. 
According to this idea, a robust MIMO-THP is proposed in 
this paper for time-varying channels based on Minimum 
Mean Square Error (MMSE) criterion where it is assumed 
that the perfect CSI is only available at the beginning of each 
transmission block, and for the next symbols, the channel 
deviations are pre-equalized blindly. Moreover, the transmit 
power of the proposed robust MIMO-THP is reduced by 
using the Tilted constellation concept. This reduction in 
transmit power in time varying channels is translated into a 
reduction of Inter Channel Interference (ICI) between MIMO 
sub-channels and in addition compensates some of the 
THP's power loss. Our simulation results show that the 
proposed robust MIMO-THP outperforms the conventional 
MIMO-THP technique over correlated channels. 
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Introduction 

It is well known that the Multiple-Input Multiple- 
Output (MIMO) systems have higher capacity than the 
traditional single-input single-output (SISO) systems; 
however, to achieve this capacity, it is required to 
cancel the interference between sub-channels at the 
receiver side. To do this, some techniques as 
Tomlinson-Harashima Precoding has been proposed 
to pre-eliminate this interference [Windpassinger 
(2004)]. In general, the ability to reduce the sub- 


channels interference depends on the available 
amount of Channel State Information (CSI) at the 
transmitter side. The common way to obtain the CSI at 
the transmitter side is to estimate them at the receiver 
and feed them back to the transmitter through a 
feedback channel. However, this method has only 
good performance in flat fading MIMO channels, and 
it is not suitable for the case of time-varying MIMO 
channel [KhaleghiBizaki et al. (2010)] in which the 
THP parameters should be designed in each (or some) 
symbol time slot. However, this method seems to be 
very complicated and tends to produce high feedback 
data overhead. Therefore, in most of the papers, the 
channel is assumed to be static during one or several 
symbol intervals, but this assumption is less valid for 
practical scenarios and may reduce the performance of 
the system [KhaleghiBizaki et al. (2010)]. 

In some works, the time-varying nature of MIMO 
channels is discussed. In Fischer et al. (2002), the 
authors proposed a combined linear precoder and 
Space-Time Block Code (STBC) method which is 
designed for rapidly time-varying channels based on 
Minimum Mean Square Error (MMSE) criterion; 
however, their method has a complicated structure 
and needs many calculation to compute the combined 
precoder and STBC. 

In Tong et al. (2011), the proposed linear precoding for 
time-varying channels based on the Conjugate 
Gradient (CG) method could reduce complexity of the 
system compared with the conventional approaches, 
while increasing the MSE. 

In Kung et al. (2004), a recursive QR approach has 
been presented for designing of linear equalizer in 
time- varying frequency selective MIMO channels. 
However, their works need many computations and 
iteration to compute precoder coefficients. 

In contrast to the previous works in which the linear 
precoder was taken into consideration, in Castro et al. 
(2005) an adaptive nonlinear precoder (i.e. the THP 
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precoder) was proposed for the frequency selective 
channel, which is optimized according to the updated 
CSI in every symbol time slot. This updating for each 
time slot needs a high volume of feedback data from 
receiver to transmitter. 

In this paper, we extend the MMSE Tomlinson- 
Harashima precoder design concepts of Khaleghi- 
Bizaki et al. (2008) to the time-varying MIMO channels 
and introduce a robust MMSE-THP precoder which is 
robust against channel variations. The proposed 
method needs the CSI only at the first of each data 
blocks and blindly pre-equalizes the channel during 
the block based on the time correlation property of 
channel. In addition. Tilted constellation method of 
Kang et al. (2009), Kim & Lee (2013), Carolin & 
Johannes (2013) is used to reduce the transmit power 
of the proposed robust THP. The power reduction 
mitigates the Inter Channel Interference (ICI) between 
sub-channels and in addition recovering some of the 
THP's power loss. Consequently, this ICI reduction 
and power loss recovering cause the BER of system to 
be improved. 

It should be noted that in the Tilted constellation of 
THP which is originally developed for SISO channels 
in Kang et al. (2009), by tilting or rotating the ordinary 
constellation by a proper angle, the transmit power is 
reduced. The appropriate angle is the one that 
minimizes the transmission power and it is chosen 
from a set of possible angles. This power reduction is 
greater at higher SNRs and shorter data block lengths. 
Since the transmitted power determines the slope of 
BER curves Kang et al. (2009), at these SNRs. a lower 
BER can be achieved with tilted-constellation THP. In 
our scheme, we extend this idea to MIMO channels 
and more over, for time-varying channels. 

In the proposed Tilted-MIMO-THP, we use different 
angles for tilting the constellation at each antenna. 
Hence, we need to find a set of proper angles, which 
minimize the total transmitted power of all antennas. 
Appropriate angles then are sent to the receiver 
through a feed-forward link. We consider MMSE- 
MIMO-THP in our proposed method and develop the 
Tilted constellation method for it in the time- varying 
channels. The performance of the proposed method is 
studied in different symbol block lengths and different 
normalized Doppler frequency shifts. By simulations, 
it is shown that the proposed method has a good 
improvement in performance with respect to the 
previous works. The complexity of our proposed 
method is also investigated with respect to some 


previous methods. The material in this correspon- 
dence was presented in part at Khaleghi et al. (2012) 
and Khaleghi Bizaki et al. (2013). 

This paper is organized as follows: section II describes 
the system model. Section III provides the detail of the 
proposed Tilted-MIMO-THP scheme. In section IV, 
computer simulations are conducted to evaluate 
performance of the proposed method. Finally, the 
conclusion is given in section V. 

System Model 

We consider a MIMO-THP broadcast communication 
system with n T transmit antennas and n R users each 
equipped with a single receive antenna. A block 
diagram of the proposed Tilted-MIMO-THP scheme is 
illustrated in Fig. 1 and briefly explained below. 



FIG. 1 PROPOSED MIMO-THP WITH TILTED CONSTELLATION 
IN DECENTRALIZED SCENARIO 


The n x dimensional vector a t represents input symbol 
and signal x t at time instant t, passing through 
feedback filter B t , and tilt matrix T which is added to 
the intended transmit vector to pre-eliminate the 
interference from previous users. Then the resultant 
signal is fed to modulo-operator to keep the resultant 
symbols within the original constellation. The 
boundary region of the constellation is related to the 
modulation constellation, which for rectangular M-ary 
QAM modulation is considered to be t = 2VM 
[Windpassinger (2004)]. Using tilted constellation, the 
signal constellation is further rotated by a set of 
appropriate angles, represented by the diagonal 
elements of tilting matrix T; to reduce the ICI and 
power loss of THP. 

The output signal of tilted matrix is then passed 
through a unitary feed-forward filter to remove the 
residual interference. Finally, the precoded signal is 
sent through the MIMO channel. In view of the fact 
that all interferences are take into account when the 
transmit signal is calculated and subtracted at the 
transmitter side. The receivers are left with some 
simple operations including power scaling (i.e. the 
diagonal matrix G t ), de-modulo operation(to keep the 
received symbols within the original constellation), 
constellation tilting with reverse angle and a single 
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user detection. According to Fig. 1, the received signal 
can be represented as: 

y t = H t x t + n t , (1) 

where X t E C n ™ , y t E C n R™ , H t E C n ** n T and 
1l t E C nRxl represent the transmitted, received, 
channel fading coefficients and noise vectors, 
respectively ( C denotes complex domain). The 
elements of the noise vector are assumed as 
independent complex Gaussian random variables with 
zero mean and variance o\, i.e. n t ~CN(0, <JnI nR ). 

In time-varying channel, it is assumed that the perfect 
CSI is known at the beginning of each data block, but 
not available during the block; meaning that it is 
desired to optimize THP at time t with the assumption 
that the outdated CSI, i.e. is available. H t _ T 

corresponding to the channel state x second earlier (i.e. 
the beginning of the block). For simplicity, we assume 
x = 1T S where 1 < N, andT s and 1 represent the symbol 
time and number of symbols, respectively. 

The outdated CSI H t _ T can be obtained from H t based 
on the correlated channel model as [Fischer et al. 
( 2002 )]: 

H t ~CN(p t H t _ T , (1 — pt)I), (2) 

where p t is common time-correlation of the i.i.d. time- 
varying MIMO channel and defined as 

p t = E{[H t ] i , j [H t _ T ]H} = R(T), (3) 

where R(x) depends on the channel time-variation 
model. 


As shown in Fig. 1, we use a titling diagonal matrix T 
in both feed-forward and feedback paths in the 
transmitter side as 


e +j01 0 0 ' 

T = 0 0 . (4) 

- 0 0 e +j 0 n T- 

It should be noted that this matrix is equal to an 
identity matrix in the conventional MIMO-THP 
(untitled MIMO-THP) systems. 


For MMSE criterion, the error is defined as the 
difference between the effective data vector v t and the 
data vector at the input of the decision modulo r t as 

e t = r t - v t = [T -1 G t H t F t - T -1 B]x t + n t , (5) 


where n t = T 1 G t n t is the equivalent noise. Using 
MMSE criteria, this error is minimized as: 


argmin E Ht | Ht _ x 

Gt,Bt,F t 


E a ,xll[T -1 G t H t F t - T-^ + n t || 2 

2 ^ r> ^ ' 


s.t. E ||x t || 2 < P T 

where P x is total transmitted power. Instead of solving 
(6), it is easier to use the orthogonallity principle 
[KhaleghiBizaki et al. (2008)]. In this case, the MMSE 


solution should satisfy: 

E Ht |H t _ T {e t rt H } = 0. (7) 

Hence 

Gt^rrx — B t cJ)^ rt , (8) 

where c|> rr , t = E Ht | Ht _ T {e t r t H } and 4> Xr , t = E[x t r t H ]. 

Using (1) the matrix 4> rr t can be obtained as: 

4*rr,t = E Ht | Ht _ x {e t r t H } = E Ht | Ht _ T {a|H t H t H + 0^1}. (9) 

Calculation of conditional expectation in (9) is not 
straightforward. We consider an approximate solution 
by considering the channel H t as H t = H t + A t where 

H t = H t + A t = Pt H t _ T +V(l - Pt)Et/ (10) 

where p t = R(x) = J 0 (2Tif D x)is the correlation coeffi- 
cient between the two instant t and t — x, where J 0 is 
the zeroth order Bessel function of the first kind 
and f D is maximum Doppler frequency, and E t is the 
circularly symmetric complex Gaussian matrix with 
i.i.d. entries (E t ~CN(0, 1)). Hence, relation (9) can be 
written as: 

<l>rr,t = E Ht |H t _>x(H t + A t )(H t + A t ) H + a 2 1} 

= a|H t H t H + C At + a 2 I ( 

where we assumed C At = (1 — p£)I. Also, to calculate 
4>^ r ,t, we can substitute (11) and (12) in (8), then we 
have: 

4hr,t = E[x t r t H ] = E Ht | Ht _ T {E a ,x[x t r t H ]} 

= E„ t | Ht _ x {E a ,x[x t x t H F t H H t H ]} 

= EH t |H t - x { C7 lEt H (Ht + A t ) H } 

= ^F t H H t H 

F t H = + C At + <;i)H t “ H , (13) 

where ^ . SinceF t is considered to be unitary 

matrix, we obtain 

s t s t H = (H t H t H + q + c At ) H t _H H t _1 (H t H t H + + c At ), (14) 

where we assumed S t = G^ 1 B t in which S t can be 
obtained by Cholesky factorization of (14). The 
feedback matrix B t , the feed-forward filter matrix F t 
and the diagonal scaling matrix G t can be expressed as 
follow [Kim & Lee (2013)]: 

G t = diag(l/ s n , ... , 1/ s NtNt ) 

B t = G t S t (15) 

F t = (p t H t _ x ) _1 (pt 2 H t _ x H t _ x H + q + (1 - Pt)l)S t _H 
If the channel is assumed to be quasistatic, i.e., p t -> 1, 
and the rotation matrix is assumed to be unitary, 
i.e. , T -> I , the above equation convert to what is 
considered as conventional MMSE-THP in this paper. 

Tilted MIMO-THP 

As mentioned before, in the general model of Fig. 1, 
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each element of the diagonal tilting matrix T 
determines the rotation angle of its corresponding 
antenna signal. 

In conventional THP technique, the main constellation 
is only extended by a modul operator on its boundary 
regions, however, in the tilted-THP, the signal 
constellation is further rotated by a set of appropriate 
angles. In order to do this, a block of symbol with 
length N is divided into n x groups (equal to the 
number of transmitting antennas ) each of length 
N/n x ; Then, every group of symbols which are 
transmitted from one of the antennas, the transmitter 
chooses an appropriate angle from a set of possible 
angles in such a way that the transmitted power is 
minimized based on some criteria. From Kang et al. 
(2009), the ordinary constellation is tilted by Q 
possible angles as 0 q , q = 1 , Q .Then, the optimal 
angle 0* is selected based on the following equation: 

q* = argmin q=1 Q ||Mod[a - ie _i0 ‘i]|| / (16) 

where a is the group of symbols in each transmitting 
antenna and i illustrates the interference sequence at 
each antenna due to the previous antennas and is 
calculated similar to the conventional THP. The 
difference between the conventional (i.e. untitled) and 
tilted THP is in rotating of the transmitted signal of the 
previous antennas by their optimal angles. 
Transmitted power and resulting optimal angle of the 
each antenna are calculated as [Kang et al. (2009)]: 

Ptilted(Q) = E [ q PJ in Q X q] 

1 ^ " 2 ( 17 > 
Xq = NzJ M0d t a_ie )0q ^ 

i=l 

Symbols a 1 transmitted from first transmit antenna are 
not affected by any interference signal, that's why we 
transmit a x without any interference, with optimal 
tilted angle 0! that results from the minimized 
transmitted power for first antenna as: 

x x = ai ej 01 . (18) 


For a 2 , which is transmitted from second antenna, the 
optimal transmit angle is achieved by computing the 
minimum transmitted power similar to the a 1 but here, 
there is interference due to the a x and hence 

x 2 = (Mod(a 2 - b^e - ' 02 ))^ 02 . (19) 

Similarly, for the kth antenna, we have 


x k = 





e j0 k (20) 


where 0 k is the optimal tilted angle for the kth antenna. 


At the receiver side, the received signal for the kth 
antenna can be shown as, 

r k = (b k x + n k ) e"j 0 k (21) 

Where b k denotes the kth row of matrix B t . By using 
(20), relation (21) can be written as. 


rk = ( x k + ^ + n k j e 


Since b kk = 1, hence after modulo operator we obtain 


a = Mod(r k ) = Mod | (a k e j0k - ^ b kj Xj + 1 


+ n k )e" 


= Mod((a k e j0k + n k )e -j0k ) 

= Mod(a k + n k e“^ 0k ) 

As shown in the tilted constellation case, the matrix B t 
is similar to the untited-MIMO-THP case; and the 
difference is only in rotation of the phase of the noise. 
Since the noise has circularly symmetric complex 
Gaussian distribution, therefore the relations and 
results for the decomposition of matrices in the tilted 
constellation case are similar to the conventional case. 

If the tilted constellation is utilized to detect the 
original symbol directly, the performance will be 
degraded due to rotation and reduction of minimum 
distance between symbols. Hence we should send the 
optimal angles to the receiver to detect the desired 
symbols. In the receiver, the symbol demodulator at 
first reversely rotates the received signal for each 
antenna and then detects symbols in a way similar to 
the conventional THP. 

Now, we discuss the complexity of the proposed 
method. It is realized that the proposed method add 
some parts to the conventional THP method. To 
compare the complexities of the proposed method and 
the conventional THP, we only consider those parts 
that are changed in our method. Any plus and minus 
operation is assumed as one calculation. For example, 
for the case of n R = n x = 4 , the complexity of the 
conventional THP is 15n at the transmitter side and 8n 
at the receiver side while for the proposed Tilted-THP 
method this complexity is 25n at the transmitter side 
and 12n at the receiver side , where n = N/n x is the 
length of the input symbol to each antenna. As it can 
be seen that the increased complexity is suitable. On 
the other hand, as shown in the next section in 
simulation results, the performance gain that we attain 
by the proposed method is more considerable with 
respect to the above complexity increasing. 
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Simulation Results 

We assume a MIMO-THP system with 4 transmit and 
4 receive antennas and 4-QAM modulation; then use 
(10) to simulate and model the channel from outdated 
CSI; as well normalized Doppler frequency f D T = 
0.001 andf D T = 0.005 and f D T = 0.01 for modeling of 
the slow, medium and fast time varying channel, 
respectively [KhaleghiBizaki et al. (2010)]. 

Fig. 2 shows the transmitted power of MMSE-THP 
system versus the number of tilted angles. In this 
figure, the SNR is assumed to be 20 dB. It can be seen 
that by increasing the number of available tilted angles, 
the transmitted power reduces until a constant as floor. 
This saturation is done due to that the matrix F is 
supposed to be unitary. Also, it can be seen that by 
decreasing the block length, the power reduction will 
be more and the power reaches faster its floor. For 
example, the transmitting power for the conventional 
MIMO-MMSE-THP is 4.91 but for the tilted 
constellation and with the similar condition, that 
decreases to 4 and 4.24 when the block length 
increases to5 and 10, respectively. 



FIG. 2 TRANSMITTED POWER VS. THE NUMBER OF TILTED 
ANGLES IN DIFFERENT BLOCK LENGTH N=5,10,20. 
(SNR=20dB, Nt=Nr=4, 4QAM). 


Fig. 3 shows BER vs. SNR for several block length and 
different normalized Doppler frequencyin theMIMO- 
MMSE-THP. It can be observed that the performance 
is improved in all block lengths compared to that in 
the conventional case, especially for high SNRs. This 
figure shows that the power reduction is greater for 
smaller blocks lengths. In addition, for the smaller 
normalized Doppler frequency, BER reduction is more. 
The result demonstrates that our proposed method 
has better performance in slow time-varying channels 
with respect to the fast time-varying channels. 



(a )f D T = 0.001 




FIG. 3 BER VS. SNR AND SEVERAL BLOCK LENGTH AND 
NORMALIZED DOPPLER FREQUENCY IN MIMO-MMSE-THP 

Conclusions 

In this paper, we proposed a robust MIMO-THP 
scheme for time-varying channels based on MMSE 
criterion, which uses outdated CSI and correlation 
characteristic of channels to calculate the precoder 
coefficients. In the proposed robust method, the 
transmit power was minimized by using the Tilted 
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constellation method. The power minimization causes 
ICI reduction between sub-channels and recovers 
some of the THP's power loss. Consequently, this ICI 
reduction and power recovering result led to BER 
improvement, especially at high SNRs. The 
transmission power has reduced further by utilizing 
more tilted angles and smaller symbol block length. In 
addition, it was demonstrated that the achieved 
performance gain by using the proposed method is 
more considerable with respect to its higher 
complexity so as to avoid changing the complexity to 
n 2 or higher order. The tilted angle sends to the 
receiver so as to be known for the receiver. 
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